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ABSTRACT 

We present new l/-band differential photometry and radial-velocity measurements of the unevolved 
1.18-day period F-|-G-type double-lined eclipsing binary VZ Cep. We determine accurate values for the 
absolute masses, radii, and effective temperatures as follows: Ma = 1.402 ± 0.015 M©, i?A = 1.534 ± 
0.012 Rq, Teff = 6690±160K for the primary, and Mb = 1.1077±0.0083 Mq, i?B = 1.042±0.039 Rg, 
Toff — 5720 ± 120 K for the secondary. A comparison with current stellar evolution models suggests 
an age of 1.4 Gyr for a metallicity near solar. The temperature difference between the stars, which 
is much better determined than the absolute values, is found to be '--^250 K larger than predicted 
by theory. If all of this discrepancy is attributed to the secondary (which would then be too cool 
compared to models), the effect would be consistent with similar differences found for other low- 
mass stars, generally believed to be associated with chromospheric activity. However, the radius of 
VZ Cep B (which unlike the primary, still has a thin convective envelope) appears normal, whereas 
in other stars affected by activity the radius is systematically larger than predicted. Thus, VZ Cep 
poses a challenge not only to standard theory but to our understanding of the discrepancies in other 
low-mass systems as well. 

Subject headings: binaries: eclipsing — stars: evolution — stars: fundamental parameters — stars: 
individual (VZ Cep) 



1. INTRODUCTION 

VZ Cephei (also known as BD +70 1199 and 
GSC 04470-01334; a = 21*^ 50"" 11!14, 5 = +71° 26' 38'.'3, 
J2000.0; V — 9.7 2) was discovered ph otographically as 
a variable star bv lGengler et al.l (|1928f l. who cla ssified it 
to be of type "Is?", a rapid irregular variable. ICannonI 
(jl934) made the first spectral type assig nment as GO. It s 
discovery as an eclipsing binary is due to lRossigerl (|1978D . 
who determined a period of 1.18336 days and showed it 
to have unequal rn inima. The system was included by 
iLacvl (|1992L l2002af ) in his photometric surveys of eclips- 
ing b inary st ars. Because of its relatively late spectral 
class. [Poppeij |l996) had it as a target in his program to 
search for late-type (F-K) eclipsing binary stars. He con- 
cluded on the basis of 4 spectrograms that both stellar 
components were likely hotter than GO. No determina- 
tions have been made of the physical properties of the 
stars, and the system has generally been neglected ex- 
cept for measurements of the times of eclipse made by a 
number of investigators since 1994. 

We began our investigation for the same reason Pop- 
per did: to test theoretical predictions of the properties 
of low-mass stars. We and other authors have previously 
found that in some of these binary systems the abso- 
lute properties are not well described by standard stellar 
evolu tion theorv ( see, e.g.. |Pot)peij|1997t IClausen et al.l 
[1991 Ribas 2Q0i iTorres et al.l 120061) . We find below 
that VZ Cep also shows some anomalies compared to 
standard models, even though its components are both 
more massive than the Sun. 

2. ECLIPSE EPHEMERIS 
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Photometric times of minimum light of VZ Cep avail- 
able from the literature are collected in Table [TJ Eclipse 
ephemerides determined by weighted least squares sepa- 
rately from the 15 primary minima and the 13 secondary 
minima gave the same period within the errors. A joint 
fit of all the measurements was then performed enforc- 
ing a common period but allowing the primary and sec- 
ondary epochs to be free parameters. Scale factors for 
the internal errors were determined by iterations sep- 
arately for the two types of measurements in order to 
achieve reduced values near unity. This solution re- 
sulted in a phase difference between the two epochs of 
/^(j) = 0.50033 ± 0.00022, not significantly different from 
0.5. For our final ephemeris we imposed a circular orbit, 
and obtained 

Min I (HJD) = 2,452,277.324478(59) + 1.183363762(84) • E. 

The uncertainties of the fitted quantities in terms of the 
least significant digit are shown in parentheses. The final 
scale factors for the published internal errors were similar 
to those in the previous fit, and were 1.38 for the primary 
and 2.42 for the secondary. We adopt this ephemeris for 
the spectroscopic and photometric analyses below. 

3. SPECTROSCOPIC OBSERVATIONS AND ORBIT 

VZ Cep was placed on the observing list at the 
Harvard-Smithsonian Center for Astrophysics (CfA) in 
2003 January, and was observed until 2007 June with an 
echelle spectrograph on the 1.5-m Tillinghast reflector at 
the F. L. Whipple Observatory (Mount Hopkins, AZ). A 
single echelle order 45 A wide centered at 5188.5 A was 
recorded with an intensified Reticon photon-counting 
diode array, at a resolving power of AA/A « 35,000. The 
strongest lines in this window arc those of the Mg I b 
triplet. A total of 39 spectra were obtained with signal- 
to-noise ratios ranging from 19 to 47 per resolution ele- 
ment of 8.5 kms^^. 
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Radial velocities were measured with the two- 
dimensional cross-cor relation technique TODCOR 
(jZucker fc Mazehlll994D . Templates for the primary and 
secondary were selected from a large library of synthetic 
spectra b ased on model atmosph e res by R. L. Kurucz ^ 
(see also iNordstrom et all Il994t iLatham et all l2002[ ). 
These calculated spectra cover a wide range of effective 
temperatures (Toff), rotational broadenings (wsini when 
seen in projection), surface gravities (logg), and metal- 
licities. Solar metallicity was assumed throughout, along 
with initial values of \ogg = 4.5 for both components. 
The temperatures and rotational velocities for the 
templates were determined by running extensive grids 
of two-dimensional cross-correlations and seeking the 
best correlation value averag ed over all e xposu res, as 
described in more detail by iTorres et all (|2002l ). The 
secondary component in VZ Cep is some 5 times fainter 
than the primary, and we were unable to determine its 
temperature independently. We therefore adopted the 
results from other estimates described below, and chose 
the nearest value in our library, which is 5750 K. Due to 
the narrow wavelength range of our spectra the derived 
temperatures are strongly correlated with the assumed 
surface gravities. The secondary log 5 presented in 
§[5] is very close to the value we assumed, but the 
primary log (7 is intermediate between 4.0 and 4.5, so 
we repeated the calculations above using the lower 
value, and interpolated. The results for the primary are 
Toff = 6690± 150 K and vsini = 57±3 kms-\ and for 
the secondary we obtain v sin i = 50 ± 10 km s~^ . Radial 
velocities were derived with template parameters near 
these values. The stability of the zero-point of the CfA 
velocity system was monitored by means of exposures of 
the dusk and dawn sky, and small run-t o-run correction s 
were applied in the manner described bv lLathamI (|1992D . 

Possible systematics in the radial velocities that may 
result from residual line blending in our narrow spectral 
window, or from shifts of the spectral lines in and out of 
this window as a function of orbital phase, were investi- 
ga ted by perform ing numerical simulations as described 
bv iTorres et all ([1997, 2000). Briefly, we generated ar- 
tificial composite spectra by adding together copies of 
the two templates with scale factors in accordance with 
the light ratio reported below, and with Doppler shifts 
for each star appropriate for each actual time of obser- 
vation, computed from a preliminary orbital solution. 
These simulated spectra were then processed with TOD- 
COR in the same manner as the real spectra, and the 
input and output velocities were compared. Experience 
has shown that the magnitude of these effects is diffi- 
cult to predict, and must be studied on a case-by-case 
basis. Corrections were determined for VZ Cep based 
on these simulations and were applied to the raw veloci- 
ties. The corrections for the primary star are small (un- 
der 1 kms~^), but for the secondary they are as large as 
11 kms~^, and as expected they vary systematically with 
orbital phase or radial velocity (see Figure [1]) . Similarly 
large corrections have been found occasionally for other 
systems using the same instrumentation (e.g.. AD Boo, 
GX Gem: [Clausen et al.l[2008l : lLacv et al.ll2008[) . The im- 
pact of these corrections is quite significant for VZ Cep: 
the minimum masses increase by 4% for the primary and 

^ Available at |http: //cf akuS . cf a. harvard ■ edu| 
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Fig. 1. — Corrections for systematics in the radial- velocity mea- 
surements for VZ Cep as a function of orbital phase and radial 
velocity (see text). Filled circles correspond to the primary and 
open circles to the secondary. 



1.9% for the secondary, and the mass ratio decreases by 
2.1%. The final velocities in the heliocentric frame, in- 
cluding the corrections for systematics, are listed in Ta- 
bleland have typical uncertainties of 1.3 kms~^ for the 
primary and 3.8 kms~^ for the fainter secondary. 

Preliminary single-lined orbital solutions using the pri- 
mary and secondary velocities separately indicated a 
zero-point difference between the two data sets (i.e., 
a difference in the systemic velocity 7), which is of- 
ten seen by many investigators in cases where there is 
a slight mismatch between the templates used for the 
cros s-correlations and the spectra of the real stars (see, 
e.g. , iPoppeil I2OOOI : iGriffin et all |2000[ ) . Numerous tests 
with other templates did not remove the offset.^ This 
most likely arises in our case because of stellar parame- 
ters (particularly the rotation) that fall in between the 
template parameters available in our library of synthetic 
spectra, which come in rather coarse steps of 10 kms~^ 
at the high rotation rates of VZ Cep. We therefore in- 
cluded this velocity offset as an additional free param- 
eter in the double-lined orbital fit, and we verified that 
when doing so the velocity semi-amplitudes (which de- 
termine the masses) are insensitive to the exact template 
parameters within reasonable limits, and are essentially 
identical to those resulting from separate single-lined so- 
lutions. Our final orbital fit is presented in Table [3l No 
indication of eccentricity was found, as expected for such 
a short period, so only a circular orbit was considered in 
the following. A graphical representation of the observa- 
tions and our best fit, along with the residuals, is shown 

* As a further test, solutions without applying the corrections 
for systematics described in the preceding paragraph gave an offset 
more than twice as large. 
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Fig. 2. — Radial- velocity measurements for VZ Cep (including 
the corrections for systematics described in the text) along with 
our spectroscopic orbital solution. Filled circles correspond to the 
primary, and the dotted line represents the center-of-mass velocity. 
Error bars are smaller than the size of the points. The O — C 
residuals are shown on an expanded scale in the bottom panels, 
where typical error bars are indicated in the upper right corner. 

in Figure [51 

The light ratio between the primary a nd secondary was 
determ ined from our spectra following I Zucker fc Mazehl 
(|1994f ). accounting for the difference in line blocking be- 
tween the primary and the much cooler secondary. After 
corrections for systematics analogous to those described 
above, we obtained Lb/La = 0.19 ± 0.02 at the mean 
wavelength of our observations. A further adjustment to 
the visual band taking into account the temperature dif- 
ference between VZ Cep A and B was determined from 
synthetic spectra integrated over the V passband and 
the spectral window of our observations, and resulted in 
{Lb/La)v = 0.22 ±0.02. 

4. PHOTOMETRIC OBSERVATIONS AND ANALYSIS 

Differential photometry of VZ Cep was obtained at 
the URSA Observatory on the University of Arkansas 
campus at Fayetteville, AR. The URSA Observatory sits 
atop Kimpel Hall and consists of a Meade f/6.3, 10-inch 
Schmidt-Cassegrain telescope with a Santa Barbara In- 
struments Group ST8EN CCD camera inside a Technical 
Innovations RoboDome, all controlled by a Macintosh G4 
computer in an adjacent control room inside the build- 
ing. The field of view is 20' x 30'. Images of VZ Cep 
{V ~ 9.7) and the two comparison stars GSC 04470- 
01497 {V « 9.9) and GSC 04470-01622 {V w 11.2), both 
of which are within 6' of the target, were taken with typ- 
ical integration times of 60 seconds through a Bessell V 
filter. With an overhead of about 30 seconds to down- 
load the images from the camera, the observing cadence 
was typically 90 seconds per image. A "virtual mea- 
suring engine" application written by Lacy was used to 
determine the brightness of the variable and comparison 




Fig. 3. — y-band photometric measurements for VZ Cep, along 
with our best constrained model described in the text. O—C resid- 
uals are shown at the bottom. 



stars, to subtract off the sky brightness, and to correct 
for differences in airmass. A total of 5473 images were 
gathered between 2001 March 5 and 2003 September 7. 
Differential magnitudes were formed between the vari- 
able star and the magnitude corresponding to the sum of 
the fluxes of the two comparison stars. These are listed 
in Table m and shown graphically in Figure [3] along with 
our modeling described below. Expanded views of the 
primary and secondary eclipse are given in Figure 0] and 
Figure [5l The typical precision of these measurements is 
about 0.007 mag, which is comparable to that expected 
from photon statistics ('~0.006 mag). The comparison 
stars are not known to be variable. The mean magnitude 
difference between the two was constant with a standard 
deviation of 0.0095 mag over 67 nights, which is what 
would be expected for individual magnitudes with a typ- 
ical error of 0.007 mag. A Lomb-Scargle periodogram 
analysis of the individual differences was performed to 
search for periodic signals in either star, but none were 
detected. 

We have previously found (|Lacv et al.ll2008f) that the 
URSA photometry is significantly improved by removal 
of small nightly zero-point variations. Thus 67 nightly 
corrections were made to the original magnitudes based 
on a preliminary photometric orbital fit. This procedure 
reduced the residual standard deviation by about 15%, a 
small but significant amount. Examination of these off- 
sets, which are typically smaller than 0.01 mag, revealed 
no detectable pattern as a function of orbital phase. Such 
a pattern might be expected, for instance, if there were 
perturbations in the light curve due to spots on either 
star (assuming synchronous rotation). We thus consider 
the nightly offsets to be instrumental in nature. 

The corrected photometry was fitted w ith the NDE 
model (|Etzemi981l : iPopper fc E"t^ |1981(). with aU ob- 
servations being assigned equal weight. In this model 
the stars are represented as biaxial ellipsoids, and de- 
spite the relatively large radius of the primary of VZ Cep 
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Fig. 4. — Enlarged view of Figure [3] showing the V-band pho- 
tometry for VZ Cep around the primary minimum. O—C residuals 
are shown at the bottom. 




Fig. 5. — Enlarged view of Figure|3]showing the V-band photom- 
etry for VZ Cep around the secondary minimum. O — C residuals 
are shown at the bottom. 

relative to the sepa ration (see b elow), its ellipticity of 
0.016, as defined by lEtzell (|1981D. is still well belo w the 
maximum tolerance of 0.04 (jPopper fc EtzellfTOSlh . and 
thus the model is expected to be adequate for this case. 
We return to this below. We used the JKTEBOP imple- 
mentation of lSouthworth et al.l ()2007f ) with a hnear limb- 
darkening: l aw, consistent with our findings (|Lacv et al.l 
120051 12008D that with the amount and precision of our 
data, non-linear limb-darkening laws do not improve the 
accuracy of the fits significantly. The following quan- 
tities were allowed to vary in this unconstrained solu- 
tion: the central surface brightness Jb of the secondary 
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Fig. 6. — Application of the constraint given by the spectro- 
scopic brightness ratio to the light curve fits of VZ Cep. Grids of 
solutions for fixed values of k are shown for several key parame- 
ters, along with the corresponding rms residual of the fit (cr). The 
spectroscopic value {Lb/ L\)v = 0.22 ± 0.02 is applied in the top 
left panel to determine fc, and all other quantities are interpolated 
to the same value. 

(cooler) star relative to the primary, the sum of relative 
radii ta + ?'b, the ratio of radii re/rA, the orbital in- 
clination i, the limb-darkening coefficients ma and uq, a 
phase offset, and the magnitude at quadrature. The fol- 
lowing quantities were kept fixed: the orbital eccentric- 
ity e = 0, the mass ratio from the spectroscopic solution 
q = Mq/Ma = 0.7900, and the gravity brightening expo- 
nents i/A — 0.25 and j/b — 0.36, set by the temperatures 
and surface gravities following Claret (1998). The uncer- 
tainties of the adjustable parameters were estimated with 
a Monte Carlo technique in which we generated 500 syn- 
thetic light curves, solved for the parameters, and calcu- 
lated the standard error of each parameter. This process 
yielded uncertainty estimates accurate to two significant 
digits, which is sufficient for our purposes. These "Un- 
constrained" results are given in Tabled Tests allowing 
for non-zero eccentricity and third light gave statistically 
insignificant values for those quantities. 

The F-band light ratio {Lb/La)v from this fit is 
consistent with our spectroscopic value from §[3l but 
formally less precise. In similar systems with partial 
eclipses, the accuracy of the parameters (more than their 
precision) is sometimes compromised because of strong 
correlations among variables and the relatively flat bot- 
tom o£_th£_^;f^surfece in the least-squares problem (see, 
e.g.. lAnderseiJll99l[ l. In such systems it is often the case 
that more accurate results are obtained by applying the 
spectroscopic light ratio as an external constraint. We 
have done this here by first computing a grid of solutions 
for a range of fixed values of k. We then interpolated in 
the smooth relation obtained between the light ratio and 
k to our value of {Lb/La)v — 0.22 ± 0.02, and derived 
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k = 0.680±0.030. Interpolations of all other quantities to 
this value of k were then carried out. This is illustrated in 
Figureinifor some of the key light-curve parameters. Note 
that a, the rms residual of the fit, changes very little for 
k between 0.65 and 0.75, demonstrating that the radius 
ratio cannot be determined accurately from photometry 
alone without an external constraint. The results for the 
light curve parameters from this constrained fit are listed 
in Table O and are adopted for further use. The uncer- 
tainties have been propagated directly from the error in 
the spectroscopic light ratio (see Figure [B]), and include 
also a contribution from the statistical uncertainties de- 
rived from a separate solution in which k was fixed at 
the best-fit value and all other parameters were left free. 

The fitted linear limb-darkening coefficients from this 
constrained solution tend to be somewhat smaller than 
predicted by theory. We find, for example, marginal dif- 
ferences of 0.1 i^l.Sa) for both stars compared to the 
calculations bv lvan Hammd ()1993l ). and more significant 
differences of 0.19 {2.5a) compared to the coefficients by 
[Claret (2000J. These diflFerences are similar in magnitude 
(and in this part icular case, o f the same sign) as those 
reported, e.g., bv lSouthworthI (j2008[ ). and may be due to 
shortcomings in the theoretical model atmospheres al- 
though observational errors cannot be ruled out. For 
further comparisons between theory and observati ons th e 
reader is referred to the recent work of iClaretl (l2008l| . 
Adopting coefficients from the tables by IClaretl ()2000r ) 
leads to values of the relative radii that are larger by 
1.1% for the primary and 1.7% for the secondary (1.6cr 
and 0.5tT, respectively). 

As a test of the reliability of the geometric parameters, 
we carried out solutions with two other light-curve mod- 
eling programs that are more sophisticated than the oii e 
we have used. One is the WINK program (|Woodlll972D . 
which adopts a better approximation to the stellar shapes 
as triaxial ellipsoids, rather than the simpler biaxial ellip- 
soids in EBOP, and includes a more detailed treatment 
of reflection effects. The version we used has been modi- 
fied and extended as described by'V^ (|1984l . ll986f l.rVa3 
([1985) ■ and Nordlund & Vaz (1990). The other program 
is the Wilson-Devinney c ode (WD: [Wilson fc DevinnevI 
[19711: lWilsonl[l979l \imi [199^ Ivan Hammel [2001 " in its 
most recent (2007) release, which uses Roche geome- 
try. Light curve solutions with these two codes were 
performed for a fixed value of fc = 0.680 (as closely as 
allowed by the different input quantities) to permit a 
direct comparison with our constrained JKTEBOP fit, 
and with the same limb-darkening law and coefficients 
as used earlier. The WINK fit delivered marginally 
smaller relative radii that differ from our JKTEBOP 
resuhs by Ava = -0.0009 and Arg = -0.0005 (i.e., 
less than 0.4%), and an inclination angle that was only 
Ai = +0°03 larger. The WD fit gave Ava = -0.0008, 
Are = -0.0005, and Ai = +0°22. These resuhs are 
thus not significantly different from those of the simpler 
model we have used, as expected from the relatively small 
ellipticity of the stars mentioned earlier. 

The individual temperatures were determined from 
the central surface brightness parameter Jb slightly ad- 
justed for limb darkening to corre spond to the disk av- 
erage (see, e. g.. iLacy et al.l [i987t ). the absolute visual 
fiux scale of [Popped ([1980D , and an estimate of the 
mean system temperature used as the initial value for 



the primary. The latter was then improved by itera- 
tion until convergence. The mean system temperature is 
based on accura t e Strom gren photometry for VZ Cep re- 
ported by [Lac^ ((2002^. Inter stellar reddening w a s esti - 
mated using the cal ibration of Pcrrv fc JohnstonI ([1982f ) 
and the method of I Crawford (1975) . which resulted in 
E{b — y) — 0.032 ± 0.007 and an intrinsic color in- 
dex of (&-?/) = 0.286 ± 0.007. The calibration by 
iHolmberg etal] ([20071 ) was then used to derive a mean 
system temperature of 6500 ± 150 K, assuming solar 
metallicity. The individual temperatures derived in this 
way are 6690 ± 160 K and 5720 ± 120 K for the pri- 
mary and secondary, respectively, which co rrespond to 
spectral types of approximately F3 and G4 ([Gravl[199^ 
p. 430). The primary Tcff is identical to our spectro- 
scopic estimate in §[3] The temperature difference based 
the light curve is of course better determined than the 
absolute values: ATeff = 970 ± 35 K. Use of a different 
color/temperature calibration for i nferring the m e an sys - 
tem temperature, such as that bv lAlonso et al.l ([1996( 1. 
yields results only ^^30 K hotter. 

5. ABSOLUTE DIMENSIONS AND PHYSICAL PROPERTIES 

The spectroscopic and photometric solutions above 
lead to the masses and radii given in Table [HI Also in- 
cluded are the predicted projected rotational velocities, 
calculated under the assumption of synchronism with the 
orbital motion. The secondary value is consistent with 
our measured v sin i from §[3l but the expected value of 
64.6 ±0.5 kms~^ for the primary seems somewhat larger 
than our spectroscopic estimate of wsinz — 57±3 kms^^. 
At face value this would indicate sub-synchronous rota- 
tion of that component, which is unexpected in a short- 
period binary such as VZ Cep. Since the primary star 
dominates the light of the system, we investigated the 
possibility that there might be a photometric signal pro- 
duced, for instance, by rotational modulation from sur- 
face features on that component. For this we examined 
the residuals from our adopted light curve solution. A 
Lomb-Scargle power spectrum did not indicate any sig- 
nificant periodicities within a factor of two of the orbital 
frequency, although the primary star is likely to be too 
hot for spots to be important (see §[7|). 

There are no measurements of the chemical abundance 
of VZ Cep. Our own spectroscopy is inadequate for this, 
and the combine d-light Strdmgren ind ices along with 
the calibration bv lHolmberg et al.[ ([2007f) indicate [Fe/H] 
= -f 0.06 ± 0.09, in which the uncertainties include pho- 
tometric errors as well as the scatter of the calibration. 
The Hipparcos catalog (iPerrvman et al.[|l997[ ) has no en- 
try for VZ Cep and no trigonometric parallax is available. 
From its radiative properties as measured here we find 
the distance to the system to be 215 ± 10 pc (similar 
distances are obtained separately for each component, 
indicating the consistency of the measured properties). 

Discrepancies described in the next section between 
our effective temperature estimates and the Toff val- 
ues predicted by stellar evolution models prompted us 
to attempt a deconvolution of the combined-light pho- 
tometry of VZ Cep, as a check on both the color ex- 
cess and the temperature difference between the com- 
po nents. We used tables of_standard Stromgren indices 
bv [Crawford[ ([1975D and [OlsenI (l984) . and synthesized 
binary stars for a range of primary indices and a fixed 
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y-band light ratio given by our spectroscopic estimate 
of {Lb/La)v — 0.22 ± 0.02. We explored a wide range 
of E{h—y) values. At each reddening we determined the 
intrinsic indices for the primary and secondary that pro- 
vide the best match t o the system v alues of 6— y, mi, ci, 
and j3 as measured bv lLacvl ()2002al ). in the sense. We 
found the best agreement for E{h — y) = 0.032, in excel- 
lent accord with our earlier determination based on the 
combined light. The measured h—y, ci, and 13 indices are 
reproduced to well within their uncertainties, and mi is 
within 1.8ct. Making use of the sam e color/temperature 
calibration bv lHolmberg et aTl (j2007f ) invoked earlier, the 
intrinsic indices for each star from this photometric de- 
convolution yield temperatures of 6670 K and 5690 K, 
once again in very good agreement with the light curve 
results. The temperature difference from this exercise is 
ATeff = 975 ± 40 K. 

6. COMPARISON WITH STELLAR EVOLUTION THEORY 

The absolute masses of VZ Cep have formal relative 
errors of 1% or better. The radius of the primary is 
similarly well determined, while that of the faint sec- 
ondary is good to about 3.7%. These values along with 
the temperatures are compared here wit h stellar evolu- 
tion models from the Yonsei-Yale series (|Yi et all 12001 
iDemargue et al.l |2004[) . In Figure [7] the measurements 
are shown in the R vs. T^s plane against evolutionary 
tracks computed for the measured masses and for solar 
metallicity [Z = 0.01812 in these models, indicated with 
solid lines). The shaded areas represent the uncertainty 
in the location of each track due to the measurement 
errors in the masses Ma and M^. While the primary 
track is in good agreement with our temperature deter- 
mination for that star, the secondary track is too hot. 
Adjustment of the chemical composition of the models to 
Z = 0.0280 (corresponding to [Fc/H] = +0.21) provides 
the fit shown with the dotted lines. This fit is marginally 
consistent with our temperature error bars in the figure, 
but the agreement is misleading since the temperature 
difference is much better determined than the error bars 
appear to indicate. The best-fit age for this metallicity 
is 1.6 Gyr, and the corresponding isochrone is indicated 
with a dashed curve. 

Figure [5] shows the measurements in the mass-radius 
diagram against the same set of models. The dashed line 
represents the same isochrone shown before, and the solid 
line is an isochrone for solar metallicity that provides 
the best fit, in this case for a slightly younger age of 
1.4 Gyr. Both are seen to represent the measurements 
equally well. 

These comparisons suggest that the models correctly 
predict the radii of the stars at the measured masses, but 
that the temperature of the secondary is underestimated 
by a signi ficant amount. Tests wit h a different series of 
models by iPietrinferni et al.l ()2004[ ) gave similar results. 

7. DISCUSSION AND CONCLUSIONS 

VZ Cep stands out among the F stars as one of the 
eclipsing binaries with the largest difference in mass be- 
tween the components {q — 0.7900), which provides ex- 
tra leverage for testing stellar evolution models. There 
are no less than five other systems with well deter- 
mined properties (BW Aqr B V1143 Cyg A BP Vul B 
V442 Cyg B, and AD Boo A: lAnderseiilU9i9lHLacv et all 




7000 6000 5000 

Effective temperature (K) 

Fig. 7 . — Stellar evolution models from the Yonsei-Yale series 
lIYi et aI.II2001l : [Demargue et al]|2004l ) compared against the mea- 
surements for VZ Cep. Solid lines show evolutionary tracks for 
the measured masses and for solar metallicity (Z = Zq), with the 
uncertainty in the location of the tracks represented by the shaded 
regions. Dotted lines correspond to mass tracks at a somewhat 
higher metallicity of Z = 0.0280 that seems to fit the observations 
better (see text). An isochrone for this metallicity and an age of 
1.6 Gyr is shown with a dashed curve. 




Mass (Mq) 

Fig. 8. — Isochro ncs from the Yonsei-Yale series l IYi et al.l 
120011 : IDemargue et aL| [2004) compared with the measurements for 
VZ Cep in the mass-radius plane. The dashed line is the same 
isochrone shown in Figure [7] (Z = 0.0280), and a solar-metallicity 
isochrone is represented by the solid curve, for a slightly younger 
age that fits the observations best. 

120031 iClausen et alll2008f) that have at least one compo- 
nent nearly identical in mass to the primary of VZ Cep 
(i.e., within 1%). However, these stars are all in very 
different evolutionary states so that their radii span a 
range of 33% and their effective temperatures differ by 
up to 360 K. They are therefore of little help in under- 
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standing the discrepancies with theory noted above for 
VZ Cep. Only one other well measured binary has one 
component with a mass similar to that of the secondary 
of VZ Cep, but that star (EK C ep B) is consi dered to be 
in the pre-main sequence stage (Popper 1987). 

Figure [7] highlights the main disagreement between the 
models and the measurements for VZ Cep, which is that 
the temperature difference predicted from theory for the 
measured masses and surface gravities is much smaller 
than all of our estimates. Solar metallicity models, which 
appear to fit the properties of the primary well, indicate 
ATeff = 710 K, and this is reduced further to 660 K 
when considering the higher metallicity oi Z = 0.0280. 
The uncertainty in these determinations is difficult to 
quantify, but a useful measure may be obtained by prop- 
agating the uncertainty in the measured masses, which 
results in ±50 K. Uncertainties from physical inputs to 
the models are unlikely to add much to this due to the dif- 
ferential nature of the comparison. In this paper we have 
made three empirical determinations of AToff , as follows: 
1) AToff = 970 ± 35 K, based on the Jb value from 
o ur light-curve a nalysis along with the visual flux scale 
of iPoppeil ()1980D and our spectroscopic brightness ratio 
(used as an external constraint); 2) ATeff = 975 ± 40 K, 
from photometric deconvolution based on the measured 
Stromgren indices and the spectroscopic brightness ratio 
(§151), along wit h the color /temperature calibrations of 
iHolmberg et"al] (|2007D : 3) AToff = 940 K, directly from 
a primary temperature estimate based on spectroscopy 
(§[3l) and an assumed temperature for the secondary sim- 
ilar to estimates for that star from the other two meth- 
ods. While these three determinations are not completely 
independent, their consistency despite the widely differ- 
ent ingredients reinforces our conclusion that the model 
AToff is at least '^250 K too small. 

Stellar evolution models have been shown previously 
to overestimate the effective temperatures of low-mass 
stars in eclipsing bi n aries b y up to ~200 K (e.g., 
[Torres fc Ribai l2002t iRibad 120 03). The study of 
V1061 Cyg by [Torres et aUmb&) suggested that the 
problem is not confined to M dwarfs, but extends up to 
masses almost as large as that of the Sun (0.93 Mq in 
the case of V1061 Cyg B). At the same time, the radii 
of these stars appear too large compared to theory, and 
both discrepancies are generally attributed to the effects 
of stellar activity in these short-period, tidally synchro- 
nized and rapidly rotating systems. 

There is little doubt that the VZ Cep system is ac- 
tive, judgin g by its strong X- ray emission as recorded 
by ROSAT IVoges et al.lfl999( ). We estimate its X-ray 
luminosity to be logTx = 30.61 ± 0.12 (where Lx is in 
cgs units), and logLx/Tboi = —3.70 ± 0.13.'"' The mass 
of VZ Cep B is slightly larger than that of the Sun, but 
it still has a thin convective envelope representing about 
1.3% of the total mass, which suggests that star may in 
fact be responsible for most of the X-ray emission given 
that the primary has no convective envelope. Another in- 
dication is given by the Rossby numbers of the stars (ra- 
tio Ro between the convective turnover time and the rota- 
tional period). For the primary we estimate log i?o > 2.1, 

^ By comparison, logLx fo r the Sun ra nges between 26.4 and 
27.7 during the activity cycle l lPeres et al.|[200ffl . and logLx/^bol 
ranges between —7.2 and —5.9. 



which according to I Hall ()1994 Fig. 6) clearly places it 
among the inactive stars. The secondary, on the other 
hand, has logi?o = —1-3. Stars in this regime tend to 
be very active and have photometric variability due to 
spots with amplitudes as large as ~0.4 mag. Detection 
of this expected variability is difficult in VZ Cep because 
of the faintness of the secondary. Nevertheless, we exam- 
ined the nightly residuals from our adopted solution near 
the primary minimum, where the contrast is more favor- 
able, and we see only occasional systematic deviations 
on one or two nights. However, similar deviations are 
seen at the secondary eclipse, and also outside of eclipse, 
which leads us to believe these are residual instrumental 
errors rather than real changes in the light level caused 
by spottedness on the secondary (see §|4]). 

If we consider the properties of the primary of VZ Cep 
to be relatively well described by theory for a metal- 
licity near solar, then the secondary shows a temper- 
ature difference compared to models in the same di- 
rection as mentioned above for other active stars (i.e., 
lower than predicted). However, we see no indication 
that its radius is larger than predicted (Figure [8]), which 
we would have expected not only from the evidence dis- 
played by other systems but also from recent theoreti- 
cal studies of the effects of chromospheric activity (e.g., 
MuUan & MacDonald 2001; Chabrier et al. 2003). 

As described in previous sections, we have carried out a 
variety of tests to explore the possibility of systematic er- 
rors in our mass, radius, or temperature determinations, 
including a careful examination of biases in our velocity 
measurements, and sanity checks of our light-curve anal- 
ysis with results from more sophisticated modeling pro- 
grams. We were unable to demonstrate any significant 
errors that would explain the discrepancies in the preced- 
ing paragraph. For example, matching the model AToff 
with the AToff measured from the light curve would re- 
quire a change in the mass ratio to q ~ 0.71, much lower 
than allowed by the spectroscopy, regardless of the choice 
of cross-correlation templates (see §[3]). Conversely, de- 
riving a smaller temperature difference from the light 
curve to match the models would require an increase in 
Jb to values that are unrealistic and would bring strong 
disagreement with the light ratio from spectroscopy. Ad- 
ditionally, this would leave the other two empirical esti- 
mates of AToff unchanged, and a discrepancy would re- 
main. As indicated earlier, we see no evidence for third 
light at a level that would make much difference. The ad- 
justments required in each of the quantities mentioned 
above, and others we experimented with, are so large 
compared to the uncertainties that a combination of ef- 
fects is unlikely to resolve the issue either. 

At the suggestion of the referee, we show in Fig- 
ure [9] a comparisoii with an alternate set of models by 
iBaraffe et al.l ()1998f ). which allows us to explore the ef- 
fect of differences in the mixing length parameter Q!ml- 
Previous studies of the radius and temperature discrep- 
ancies for active low-mass stars have indicated that a 
value of ttML lower than appropriate for the Sun, repre- 
senting a reduced overall convective efficiency, provides a 
much better match to the observations of these objects. 
Consequently, we show a solar-metallicity evolutionary 
track for a solar-like mixing length parameter for the ra- 
diative primary (aML = 1-9 in these models), and tracks 
for three values of the mixing length parameter for the 
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Fig. 9. — Radius and effective temperature of VZ Cep compared 
agains t evolutionary tracks for solar metallicity by IBaraffe et al.l 
II1998I) . for the exact masses we measure (solid lines). A single 
track is shown for the primary star, for a mixing length parameter 
c«ML = 1-9 appropriate for the Sun. The dashed line represents 
the same solar-metallicity Yonsei-Yale track shown in Figure [7] 
Three Baraffc et al. ( 199g) models are shown for the secondary, for 
different values of QmLi ^ labeled. For reference, the triangles on 
the solid curves correspond to an age of 1.4 Gyr, which was found 
in Figure|8]to provide a good match in the mass-radius plane using 
solar-metallicity models from the Yonsei-Yale series. 

secondary star (q;ml = 1-9, 1.5, and 1.0), which, as men- 
tioned earlier, we believe to be the more active member of 
the system. For reference, triangles on each track mark 
the age of 1.4 Gyr, which we found in Figure[7]to provide 
the best match for Z = Zp, u sing the Yonsei-Yale mod- 
els. The IBaraffe et aD (jl998[ ) model for the primary is 
seen to be similar to the corresponding solar-metallicity 
Yonsei-Yale model (dashed line, reproduced from Fig- 
ure [7]) . A reduction of the mixing length parameter for 
the secondary star leads to the expected systematic de- 
crease in effective temperature, and an increase in ra- 
dius. A secondary model with aML between 1.0 and 1.5, 
when paired with the standard aML — 1-9 model for the 
primary, would appear to give approximately the cor- 
rect temperature difference for the system. However, the 
measured radius of VZ Cep B is somewhat smaller than 
predicted, in agreement with our earlier conclusion that 
this star appears normal in size (compared to standard 
models). 
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the radius. VZ Cep thus presents a challenge to theory. 
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whether Ca II emission is present. 
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TABLE 1 

Published measurements of the times of eclipse for VZ Cep. 



HJD 




Uncertainty 


(O-C) 




(2,400,000+) 


Typc^ 


(days) 


(days) 


Source 


dQ'^f\7 4991 




0, 


.0003 




\ 


•^1 fins 79^^70 




0, 


,0003 




2 
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2 


0, 


.0005 


_n onoi 2 


3 


'^9^4zL 7Q41 


2 


0, 


.0005 


4-0 ooo'^fi 


3 


'^9n'^1 8Q4in 


2 


0, 


,0005 


4-0 0001 8 


3 


.'19054 8591 5 




0, 


.00011 


4-0 00008 


3 


5907*^ 78570 




0, 


,0002 


— 0001 Q 


3 


5907fi 74440 


2 


0, 


,0003 


— 0001 (\ 


3 


59n7Q 7(^^00 




0, 


,0003 


-UO 0009Q 


3 


59080 88fi04 




0, 


,00019 


— 0000'^ 


3 


590QS Q0970 




0, 


,0005 


— 000*^8 


3 


591 08 fiQfi'^0 


2 


0, 


,0003 


-UO 000Q9 


3 


591 1 1 fi5970 




0, 


,0006 


— OOOSS 


3 


591 1 9 8'^70Q 




0, 


,00014 


4-0 0001 Q 


3 


591 1 4 fil 500 


2 


0, 


,0006 


-UO 00980 


3 


521 54 S447n 


2 


0, 


,0004 


—0 001 87 


3 


52159 58000 


2 


0, 


,0010 


—0 00002 


3 


52166.67970 


2 


0, 


,0003 


-0.00051 


3 


52179.69710 


2 


0, 


,0003 


-0.00011 


3 


52233.54070 


1 


0, 
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+0.00070 


3 


52277.32429 


1 


0, 


,00007 
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4 


52463.70530 


2 


0, 
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+0.00079 


5 


52464.88810 


2 


0, 


,0003 


+0.00022 


5 


52482.63870 


2 


0, 
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+0.00037 


5 


52518.73064 


1 


0, 


,00019 


-0.00003 


5 


53366.01950 


1 
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+0.00037 


6 


53658.30900 


1 


0, 


,0006 


-0.00098 


7 


54009.76910 


1 


0, 


,0001 


+0.00008 


8 



Note. — Refe renc es: 1. lAgerer &: Huebscheii II1995I) : 2. INelsonI 1 120011) : 3. ILacv et all I I2002I) : 4. ISarounova fc Wolj I l2005f) : 5. ILacvl ll2002bD : 6. 
IKim et al.l ll2006f) : 7. IDiethlemH2006D : 8. INelsonI IfoOTI) . 
^ Type: 1 — primary eclipse; 2 — secondary eclipse. 
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TABLE 2 

Radial velocity measurements of VZ Cep. 



HJD 


RVa 




(U-C)A 


(C-CjB 




(2,400,000+) 


(kms ^) 


ti — 1 \ 
(kms ^) 


(kms ^) 
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+1.01 


0.692 




06 ^6 


— 150.44 


—0.27 


+4.89 


0.720 




— 109.05 


+117.57 


—0.46 


+0.23 


0.344 


53336.6417 


-116.21 


+124.74 


-0.51 


-1.60 


0.175 


53339.6764 


+108.41 


-157.55 


-0.30 


+0.16 


0.7.39 


53630.7663 


+106.36 


-156.85 


-1.07 


-0.75 


0.724 


53636.7401 


+108.74 


-160.82 


+0.93 


-4.24 


0.773 


53690.6534 


-112.94 


+127.71 


+0.49 


+4.24 


0.332 


53691.7072 


-127.60 


+134.58 


-0.62 


-6.04 


0.222 


54042.6132 


+111.73 


-160.65 


+2.80 


-2.66 


0.755 


54043.7068 


+98.06 


-145.42 


+0.68 


-2.05 


0.679 


54074.5790 


+106.66 


-154.30 


-1.60 


+2.85 


0.768 


54103.5713 


-127.62 


+139.12 


+0.44 


-2.87 


0.268 


54279.9427 


-121.65 


+138.43 


-1.19 


+6.07 


0.310 


54282.8950 


+ 105.11 


-148.54 


+3.09 


+0.71 


0.805 



Note. — These velocities include corrections for systcmatics (see text). 
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Torres & Lacy 



TABLE 3 

Spectroscopic orbital solution for VZ Cep. 



Parameter Value 



Adjusted quantities 

P (days)=>- 1.18336377 

Ti (HJD-2,400,000)'' 52,277.32446 

Ka (kms-i) 118.88 ± 0.22 

Kb (kms^i) 150.48 ± 0.67 

7 (kms-i) -9.90 ± 0.21 

AW (kms-i)'' -2.31 ± 0.65 

Derived quantities 

MAsin3j(M0) 1.339 ± 0.013 

Mb sin3 i (Mq) 1.0577 ± 0.0064 

q = Mb/Ma 0.7900 ± 0.0038 

OA sin « {10^ km) 1.9345 ± 0.0036 

ae sini (lO** km) 2.4487 ± 0.0109 

asini (Rq) 6.298 ± 0.016 

Other quantities pertaining to the fit 

A^obB 39 

Time span (days) 1621.3 

o-A (kms-i) 1.27 

o-B (kms^i) 3.82 



^ Ephemeris adopted from §[2]^ Velocity offset in the sense (primary — secondary) (see text). 

TABLE 4 

Differential F-band magnitudes of VZ Cep. 



HJD-2,440,000 AV Orbital phase 



51973.95823 0.129 0.64077 

51973.95926 0.115 0.64164 

51973.96028 0.123 0.64250 

51973.96130 0.126 0.64336 

51973.96233 0.126 0.64423 



Note. — Table^is available in its entirety in the electronic edition of the Astronomical Journal. A portion is shown here for guidance regarding 
its form and contents. 



TABLE 5 

Photometric orbital solutions for VZ Cep. 



Parameter Unconstrained fit Constrained fit 



Jb 0.495 ± 0.010 0.4920 ± 0.0013 

k = tb/t-a 0.717 ± 0.041 0.680 ± 0.030 

rA + rB 0.4077 ± 0.0043 0.4028 ± 0.0077 

ta 0.2374 ± 0.0033 0.2398 ± 0.0017 

rB 0.1703 ± 0.0073 0.1630 ± 0.0061 

UA 0.499 ± 0.075 0.420 ± 0.076 

ub 0.581 ± 0.059 0.500 ± 0.076 

i (deg) 79.47 ± 0.44 79.97 ± 0.45 

La{V)'' 0.802 ± 0.022 0.820 ± 0.014 

{Lb/La)v 0.246 dz 0.033 0.220 ± 0.020*' 

cry (mmag) 7.4400 7.4365 

A^obs 5473 5473 



^ Fractional luminosity of the primary. Adopted as a constraint from spectroscopy (sec text). 



VZ Cep 
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TABLE 6 
Physical properties of VZ Cep. 



Parameter Primary Secondary 



Mass (Mq) 1.402 ± 0.015 1.1077 ± 0.0083 

Radius (R©) 1.534 ± 0.012 1.042 ± 0.039 

logs (cgs) 4.2130 ± 0.0080 4.446 ± 0.033 

Temperature (K) . . . . 6670 ± 160 5720 ± 120 

logL (Lq) 0.634 ± 0.041 0.026 ± 0.050 

i;sini (kms-i)a 57 it 3 50 ± 10 

i)syncsini (kms-i)'' . 64.6 ± 0.5 43.9 ± 1.6 

a (R©) 6.396 ± 0.019 

Distance (pc) 215 ± 10 

Mboi (mag) 3.18 ± 0.10 4.68 ± 0.12 

My (mag) 3.16 ± 0.11 4.77 ± 0.12 



^ Value measured spectroscopically. Value predicted assuming synchronous rotation. 



